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Diamagnetism above in underdoped Bi2.2Sri 8Ca2Cu30io+(5 

Tetsuya Iye^*''', Takenori Nagatochi^^, Ryusuke Ikeda^, and Azusa Matsuda^ 

^Department of Physics, School of Science and Engineering, Waseda University, Shinjuku, Tokyo 169-8555 
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Single crystals of Bi2+a;Sr2-xCa2Cu30io+5(Bi2223) with x = 0.2 were grown by a traveling 
solvent floating zone method in order to investigate the superconducting properties of highly 
underdoped Bi2223. Grown crystals were characterized by X-ray diffraction, DC susceptibility 
and resistivity measurements, confirming Bi2223 to be the main phase. The crystals were an- 
nealed under various oxygen partial pressures to adjust their carrier densities from optimally 
doped to highly underdoped. The fluctuation diamagnetic component above the supercon- 
ducting transition temperature Tc extracted from the anisotropic normal state susceptibilities 
Xah{T) (H 1. c) and Xc{T) (H \\ c) was found to increase with underdoping, suggesting a 
decrease in the superconducting dimensionality and/or increase in the fluctuating vortex liquid 
region. 
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1. Introduction 

In multilayered cuprate superconductors, which have 
more than three Cu02 layers in their unit cells, the su- 
perconducting transition temperature Tc takes its high- 
est value when the number of Cu02 layers n is 3.^'^^ The 
Tc reduction for increasing n can be understood in terms 
. of inequivalent doping among the layers. However, the 
mechanism which enhances the Tc for the first few layers 
. is still unknown and may be related to the mechanism of 

■ high-Tc superconductivity. 

Nuclear magnetic resonance (NMR), and angle re- 
' solved photoemission spectroscopy (ARPES) have shown 
I evidence of the inequivalent doping among layers. The in- 

■ ner Cu02 planes (IPs) with square oxygen coordination 
. have lower hole doping levels than the outer Cu02 planes 
' (OPs) with pyramidal oxygen coordination. ^^^^ In the 

NMR experiments, the remarkable coexistence of antifer- 
romagnetic (AFM) IPs with superconducting (SC) OPs 
has been reported in the underdoped regime. ^^^^^ It has 
been shown theoretically that high Tc superconductivity 
is possible in this AFM and SC coexisting state. In 
that case, a finite Josephson coupling between SC OPs 
through the AFM IP provides bulk superconductivity. 
' Thus, physical properties of underdoped region of multi- 
layered superconductor is intriguing. 

In the overdoped to slightly underdoped region, sev- 
eral anisotropic properties of a trilayered superconduc- 
tor Bi2Sr2Ca2Cu30io+5 (Bi2223) have been measured 
and a sizable change in anisotropy was reported. The 
anisotropy was lower than that of double layered super- 
conductor Bi2Sr2CaCu20io+(5 (Bi2212), indicating that 
the inter layer coupling of the three Cu02 layers exists. 
Bi2223 (Tc,max = 110 K) is one of the few multilayered 
cuprates of which we can obtain large size single crys- 
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tals in order of millimeters ^^^^^ and is the simplest sys- 
tem among multilayered cuprates. In this paper, we first 
report on a growth of Bi2223 crystals suitable for un- 
derdoping. Then, the doping dependence of fluctuation 
diamagnetism in Bi2223 is extracted from magnetiza- 
tion measurements. The extracted diamagnetism above 
Tc is expected to be qualitatively the same as the dia- 
magnetism which Li et al observed through Nernst and 
torque magnetization experiment. 

2. Experimental 

Two techniques were used to reduce the doping level 
and drive Bi2+a^Sr2-a^Ca2Cu30io+5 into the highly un- 
derdoped regime. The first is to anneal samples under 
various oxygen partial pressures. Oxygen annealing con- 
trols the excess oxygen content (5, which supplies hole 
carriers to the Cu02 plane. The hole doping level can be 
decreased by lowering the oxygen partial pressure. The 
second method is a substitution of the trivalent cation 
g^3+ ^Yie divalent cation Sr^+. The hole density is 
expected to decrease with Bi substitution if the oxygen 
content is unchanged. In Bi2212, a hole density reduction 
was confirmed through Tc reduction due to this substi- 
tution. Bi2+a;Sr2-a;Ca2Cu30io+(5 single crystals with 
X = 0.1 can be grown by a traveling solvent floating zone 
(TSFZ) method. ^^"^''^ However, no successful single crys- 
tal growth of Bi2223 with x > 0.1 has been reported, 
because the growth is difficult even for the intrinsically 
substituted x = 0.1 sample. Here, using almost the same 
technique as for x = 0.1, we succeeded in growing crys- 
tals of Bi2223 with x = 0.2. 

Powders of Bi203, SrCOs, CaCOs, and CuO (all of 
99.9 % or higher purity) were mixed in the desired cation 
ratio Bi : Sr : Ca : Cu = 2.2 : 1.8 : 2 : 3, ground in an 
agate mortar, and then calcined at 760 °C for 12 hours. 
The calcined powders were well reground and again cal- 
cined at 780 °C for 12 hours. After that, the powders 
were hydrostatically pressed into a cylindrical rod under 
40 MPa and then sintered at 860 °C for 50 hours. The 
sintered rod was cut into the long (~ 7 cm) and the short 
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Fig. 1. (Color online) (a)Tlie grown crystal rod of 47 mm in 
length and typical diameter 5.5 mm. The linear scale is in cen- 
timeters. (b)One of the Bi2223 single crystals cleaved from the 
rod shown in (a). 



2 cm) part. The long one was hung on the upper shaft 
of an infrared radiation furnace equipped with two ehip- 
soidal mirrors (NEC-SCE-EDH). The short one was held 
on the lower shaft. The long one was pre- melted using 
the short one as a base for the crystal growth at a rate of 
25 mm/h in air. The pre- melted rod was cut into the pre- 
melted part (~ 6 cm) and the base (seed) part (~ 2 cm). 
The pre-melted and the seed rod were held on the upper 
and the lower shafts respectively, and both shafts were 
counter-rotated at a rate of 11 rpm. A crystal rod was 
grown using a very slow rate of 0.04 mm/h for the pre- 
melted rod. In order to realize as fast a crystal growth 
rate as possible, 300 W halogen lamps were used as light 
sources, which gave a steep temperature gradient near 
the soHd-liquid interface. It also prevents the swells of 
soHd-liquid interfaces on both the pre-melted and grown 
crystal sides which disturb static growth. Under these 
conditions, a crystal rod of 47 mm length and 5.5 mm in 
diameter as shown in Fig. 1(a) was grown. 

From the grown crystal rod, one crystal with maximum 
dimensions up to 2mm x 3mm x 0.1mm was cleaved and 
is shown in Fig. 1(b). Each part of the as-grown rod was 
evaluated by X-ray diffraction (XRD), magnetizaton and 
resistivity measurements. XRD measurements were done 
on a diffractometer (Rigaku RINT-TTRIE) using a con- 
ventional 6 — 26 method to check phase purity and deter- 
mine c-axis length. Tc was determined from magnetiza- 
tion measurements under a magnetic field of 4 Oe using a 
SQUID magnetometer (Quantum Design MPMS-7). Re- 
sistivity was measured by a standard four-probe method. 
After finishing the above sample characterizations, the 
highest quality crystal among them, whose weight was 
8.6 mg, was selected. After each annealing condition; (1) 
O2 : 760 torr, 400 °C, 72 hours, (2) O2 : 760 torr, 600 
°C, 24 hours, (3) O2 : 7.6 torr, 600 °C, 24 hours, (4) O2 
: 7.6 X 10-2 torr, 600 °C, 24 hours, (5) O2 : 7.6 x 10"^ 
torr, 600 °C, 12 hours, Tc and c-axis length were deter- 
mined, then anisotropic normal state susceptibilities Xab 
with H ± c and Xc with H \\ c were measured under a 
magnetic field of 50 kOe using the MPMS. For each an- 
nealing process, 1 hour was taken to raise temperature to 
the target value, while oxygen pressure was kept at the 
target value. At the end of each annealing step, the sam- 
ple was quenched to room temperature in 20 minutes in 
the same atmosphere in order to prevent a further oxy- 
gen uptake or loss. SCF diamagnetic components were 
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Fig. 2. (Color online) (a)X-ray diffraction peaks of the as-grown 
sample with x = 0.2. (b) Temperature dependence of DC sus- 
ceptibility for the same sample with magnetic field 4 Oe parallel 
to 0,6- axis, (c) Temperature dependence of in-plane resistivity for 
another as-grown sample. 



extracted from the normal state susceptibility data for 
each annealing condition. 

3. Results and Discussion 

Figure 2(a) shows XRD peaks in the highest quality 
c-axis oriented as-grown crystal. Only reflections from 
Bi2223 can be seen. Figure 2(b) shows the temperature 
dependence of the magnetic susceptibility in a magnetic 
field of 4 Oe parallel to the a6-plane. Clear Meissner dia- 
magnetism was observed below about 110 K. Throughout 
this report, TcS are defined as the temperatures at which 
hnear extrapolations of the steepest part of shielding sig- 
nals cross zero susceptibility. Using this definition, Tc of 
the main phase was estimated to be 110 K. A step-like 
feature at around 80 K which may attributed to Bi2212 
impurity phase was also observed in the shielding signal. 
Since XRD did not show any traces of bulk Bi2212 phase, 
it might be due to Bi2212 intergrowth, existing even in 
the highest quality crystal. By comparing the magnitude 
of the susceptibility at 80 K with the magnitude at 10 K, 
the relative volume fraction of Bi2223 phase to that of 
overall superconducting phase was estimated to be more 
than 60 %. The temperature dependence of the in-plane 
resistivity of a crystal cleaved from the same part of the 
rod was shown in Fig. 2. It showed typical T-linear be- 
havior above Tc and zero resistivity below 108.5 K. 

Next, let us discuss the precise determination of the 
c-axis lengths in the Bi2223 crystal for x = 0.2. The 
c-axis length free from systematic measurement errors 
is usually estimated by an extraporation method using a 
Nelson-Riley (NR) function (cos^ 6/ sinl9 + cos2 6/6)/2. If 
c-axis lengths obtained from (OOL) peaks i.e. dL, where 
d = A/2sin6>L, L is the Miller index of the (001) direc- 
tion are plotted against the NR function and the plots 
are fit with a straight line, an accurate c-axis length can 
be obtaind by extrapolating the fitting line to 6 = 7r/2. 
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Fig. 3. (Color online) Comparison of the relationships between 
Tc and c-axis length for x = 0.2 (filled circle) and x = 0.1 (open 
circle). The attached numbers correspond to the annealing con- 
ditions (1) to (5) for both x = 0.1 and 0.2 data. Data for a: = 0.1 
were taken from ref.-^^^ 



However, the plots for our crystal with x = 0.2 did not 
obey this function. It is known that if there is intergrowth 
of a second phase in a layered structure, dL can oscil- 
late as a function of L with a period c/Ac, where Ac 
is the difference between the c-axis lattice length of the 
main phase and that of the intergrown phase. How- 
ever, due to the strong and somewhat diverging oscilla- 
tion in the plots of our sample, it was difficult to execute a 
reliable extrapolation. Here, we computed c-axis lengths 
for X = 0.2 by averaging dL weighted by tan^/,, assum- 
ing the uncertainty of dL to be proportional to cot Ol- 
c = T^LdL tdLiiOL /^L tsLii 6 L' This corresponds to a ^ av- 
erage of dL with the weight of sin 6. We confirmed that 
this method yields nearly equivalent or slightly overesti- 
mated c-axis lengths for x = 0.1 samples if we use the 
diffraction data in ref.^^^ 

In Fig. 3, the relationship between Tc and c-axis length 
for the sample with x = 0.2 is shown by filled circles. For 
comparison, x = 0.1 data from ref.^^^ are also plotted 
with open circles. In cuprate superconductors, there is 
a general tendency that the c-axis length decreases with 
increasing hole doping due to stronger Coulomb attrac- 
tion between the negative oxygen ions and the in-plane 
hole carriers. In keeping with this tendency, our x = 0.2 
crystal shows a qualitatively similar Tc - c-axis length 
relationship to that of the x = 0.1 sample. However, sev- 
eral differences are evident. First, our sample has shorter 
c-axis lengths for the same annealing conditions except 
for the case of the condition (1), the strongest oxidation 
condition. From this, we can confirm that Bi substitu- 
tion with X > 0.1 is achieved. Second, the Tc saturation 
observed on the overdoped side for x = 0.1 seems not 
to be realized even for the same oxidation condition. For 
the X = 0.2 sample, the c-axis contracted little from the 
oxidation conditions (2) to (1), while it contracted sub- 
stantially in the x = 0.1 sample. This indicates that ox- 
idation becomes more difficult on the overdoped side in 
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Fig. 4. (Color online) Temperature and doping dependences of 
anisotropic normal state susceptibilities (a)xab(T) (H 1. c) and 
(b)xc(T) (if II c) under a magnetic field H = 50 kOe. 



the X = 0.2 sample, possibly due to the Bi substitution. 
The effect of Bi substitution in Bi2212 was studied in 
detail in ref.^^^ and summarized as follows: (1) the re- 
duction of doping level is proportional to the increase in 
Bi content, (2) the maximum superconducting transition 
temperature Tc,max is reduced, (3) there is an increase in 
excess oxygen. In Bi2212, the excess negative charges in- 
troduced to the plane by substituted Bi ions are partially 
compensated by the increase in excess oxygen and only 
20 % of substituted Bi contributes to the reduction in 
doping level. Since x = 0.2 and 0.1 samples of Bi2223 
show comparable TcS for the same oxidation conditions, 
excess Bi is thought to be fully compensated by the ad- 
ditionally absorbed oxygen. This gives a further lattice 
contraction along the c-axis as can consistently be seen 
in Bi2212.^^^ Moreover, the increase in the oxygen up- 
take due to substituted Bi makes further oxidation dif- 
ficult on the overdoped side. In addition, unexpectedly, 
Tc,max of the X = 0.2 sample is essentially same as that 
of X = 0.1 sample, while in Bi2212 about an 8 K reduc- 
tion in Tc,max was reported. This is a novel feature of 
the multilayered structure, where the bulk Tc is believed 
to be determined by the highest one among the intrinsic 
TcS of the nonuniformly doped individual Cu02 layers. 
On the overdoped side of multilayered superconductors, 
Tc is thought to be determined by the optimally doped 
jp 26,28) Disorder introduced by Bi^+ ions on the Sr site 
could only affect the OP which does not directly gov- 
ern the Tc of the sample. Thus, in this x = 0.2 case, 
the disorder might not affect the IP and Tc,max was kept 
unchanged. The invariable Tc,max indicates that the in- 
teraction between IP and BiO — SrO layer is very weak 
in optimally or slightly overdoped region. 

In Fig. 4, temperature dependences of the anisotropic 
normal state susceptibilities Xab and Xc at various doping 
levels are shown. Here, we assume that these anisotropic 
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Fig. 5. (Color online) Xab — Xc plots with an implicit parameter 
temperature T for each doping level. 



susceptibilities are composed of five components as fol- 
lows: 



vv 



XaiT) =X''- 



(1) 



where x^^^ is the isotropic Larmor diamagnetic sus- 
ceptibility, is the anisotropic Van Vleck param- 
agnetic susceptibility, Xop"a(^) XiPa(^) 
anisotropic spin susceptibilities in the OP and IP, re- 
spectively, and x™(T) is the anisotropic diamagnetic 
susceptibility due to fluctuation diamagnetism above Tc. 
Both x^^^ ciiid are expected to be small and tem- 
perature and doping independent over the experimental 
temperature range. A possible contribution to Xa from 
Landau diamagnetism is absorbed into the spin suscepti- 
bilities because it has the same temperature dependence 
as x^^^"(^)- The overall magnitudes of Xab and Xc mono- 
tonically decrease with decreasing doping level. This is 
derived from a reduction in spin susceptibility due to 
the decrease in electronic density of states (DOS) near 
the Fermi level, as well as the opening of a pseudogap 
in the under doped regime. The absolute values and the 
temperature dependences of Xab for as-grown, anneal(l), 
anneal(3), and anneal(5) in this work are comparable to 
those of Xab with x = 0.1 for corresponding annealing 
condition in ref.^^^ 

Now let us discuss the temperature dependence of Xab 
and Xc at each doping level. The gradual decrease in the 
susceptibilities can be seen from well above Tc. This de- 
crease should be caused by the opening of a pseudogap. 
However, we did not determine the characteristic tem- 
perature T* of the opening of the pseudogap since T* is 
expected to become very high in the underdoped regime 
as in the case of Bi2212. When temperature is further 
lowered to the vicinity of Tc, x™ becomes significant. 

The fluctuation diamagnetic component can be ex- 
tracted from the original data if the background nor- 
mal state susceptibility is known. However, it is ap- 



Fig. 6. (Color online) T — Tc dependence of for each doping 

level, where TcS are determined by magnetization measurement 
applying the magnetic field 4 Oe for each. The inset shows tem- 
perature and doping dependences of log|(5x^^|- Each solid line is 
drawn by fitting the data from TpD down to the temperature 15 
K lower than TpD • The absolute value of the slope of the fitting 
line corresponding to the exponent a increases with increasing 
doping level. 



parently temperature dependent in this case and there 
is no reliable way to estimate it. Here, we used the 
strong anisotropy in the fluctuation diamagnetism to 
extract it. First, Xc(T) was plotted as a function of 
XabiT) with temperature T as an implicit parameter, 
following the method described in ref.^^'^^^ The plots 
are shown in Fig. 5 for each doping level. In the re- 
gion far above Tc, they showed linear relations with al- 
most the same slope. This indicates that the tempera- 
ture dependent spin susceptibility has a doping and in- 
dependent anisotropy ratio. The deviation from the lin- 
ear background is attributed to x™ which has a differ- 
ent anisotropy ratio from that of the spin susceptibility. 
Based on the simple Gaussian superconducting fluctua- 
tion theory,^^^ anisotropy in x™ comes from anisotropy 
in the effective shape of fluctuating domain whose di- 
mension is proportional to (^ab x iab x ic-, where and 
are the in-plane and out-of-plane coherence lengths 
respectively, x™ can be estimated as the vertical devi- 
ation from the linear background assuming Xop"q;(^) ~ 



19a 



OP 



where 7 is a constant and ga is a gyromagnetic ratio of 
the conduction electrons, and that g^^ and g^^ are tem- 
perature and doping independent. Using eq.(l) and the 
above assumptions, the fluctuation diamagnetic compo- 
nent Sx^^ is calculated to be 



Xab{T) - xo 



X™(T) 



dab / 



1 xlFiT) 



(2) 
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where and Xo ^ {1 - (5c"V5aT)'}x'" + 

~ {9^^ / 9ah)'^Xah Correspond to the slope and in- 
tercept of the Hnear part of the Xah — Xc plot respectively. 
Since thought to be much smaller than Xc ^ and 

{9?^ /9ahf is less than 1 (see Fig.5), 6x^^ is nearly 
equal to X™- Thus, x™ can be estimated under the 
above assumptions. 

T — Tc dependences of Sx^^ for various doping lev- 
els are shown in Fig. 6. Here, we used T^s as the super- 
conducting transition temperature determined by mag- 
netization measurements under a magnetic field of 4 Oe. 
(5x^^s increase divergently toward Tc. We estimated the 
characteristic temprature TpD at which ^x^^ becomes 2 
% of the full value at Tc for anneal(5). Obtained Tfds 
were 134 K, 132 K, 130 K, and 116 K for samples that 
were measured as-grown and annealing under the condi- 
tions (1), (3), and (5) respectively, and they are plotted 
against 5c (the variation in the c-axis length from that 
of the as-grown sample) together with corresponding T^s 
110 K, 108 K, 99 K, and 77 K respectively in Fig.7. For 
each doping level, TpD was about 30 K higher than Tc, 
and these results are similar to the experimental results 
of Bi2212 reported by Wang et al^'^^ demonstrating that 
Sx^^ is a fluctuation diamagnetic component. Sx^^ also 
increases with decreasing doping. This behavior can be 
interpretted as that the difference between Tc and TpD 
increases with decreasing doping, x™ is expressed as 
functions of doping level p and temperature T as follows: 



■eabiP,T){r\siP,T) (3) 



Vip,T) 

where V ~ tt^^^ • mdix{£^c{p^T)^d) is the coherence vol- 
ume, d is the distance between Cu-0 blocks, and (r'^)eff is 
the mean-square radius of the fluctuating domain, which 
is roughly represented as {r'^)eff ^ £,ab for || c in this 
case. Consequently, x™ is calculated to be 



X™(p,T)«fcBT 



(4) 



Sx^^ is plotted as a function of ln(T — Tc) as can 
be seen in the inset of Fig. 6. Sx^^ exhibits a divergent 
temperature dependence (T — Tc)" with a doping depen- 
dent exponent a varing from -1.4 to -2.0. In addition, the 
exponents a determined by our experiments are rather 
lower than those expected in a classical 2D superconduc- 
tor {a = -l).3i) The same analysis on Bi2212 gave an 
essentially doping independent exponent a ~ —2.3,^^^ 
lower than that for any doping of Bi2223. Recent Nernst 
experiments for cuprates indicates that unbound vortex- 
antivortex pairs exist well above Tc and they disorder 
the long-range phase coherence. In such vortex liquid 
state, temperature dependence of the correlation length 
^ is expressed as 



C ^ exp 



B 



1/2- 



(5) 



where 5 is a constant. Since the exponential temper- 
ature dependence of is stronger than any other power 
law temperature dependence of the form of (T — Tc)", 
the nominal exponent will be lifted up if includes 
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Fig. 7. (Color online) The superconducting temperature TcS and 
the onset of fluctuation diamagnetism temperature Tp^s are 
shown with fllled triangles and open inverted triangles respec- 
tively. Tfds show the dome shaped doping dependence as well 
as TcS. 6c is the variation in the c-axis length from that of the as- 
grown sample. Since the TpD curve highly resembles the Tonset 
curve in the phase diagram of Bi2212 reported by Wang et al.,^^) 
the same physical property might have been detected in our DC 
susceptibility measurements and their Nernst experiments. 



some component of The disagreement between the 
observed exponents of Bi2223 or Bi2212 and the ex- 
pected one in simple Gaussian fluctuation theory im- 
plies that both Bi2223 and Bi2212 includes vortex liq- 
uid region above Tc. The decrease in |a| of Bi2223 on 
doping suggests decrease of the vortex liquid domain in 
the sample and/or that the system becomes even less 
two-dimensional. This behavior is associated with a more 
rapid weakening of the interlayer coupling of Bi2223 on 
underdoping than that of Bi2212. The deviation of the 
plots from fitting line near Tc may result from the fact 
that the temperature dependence of Xc ^ becomes weak 
since max(^c7<^) changed from d to due to divergence 
of on cooling down toward Tc. 

One possible origin of such modification is the exis- 
tence of the IP in Bi2223. In this system, nonuniformly 
doped Cu02 layers were confirmed by ARPES measure- 
ments, which estimated the doping level of the IP to be 
7 %, while that of OP was 23 % in optimally doped sam- 
ples. According to the ARPES result, slight underdop- 
ing can easily drive an IP insulating or possibly into an 
AF ordered state. The insulating IP disrupts supercon- 
ducting coherence among the three layers (OP-IP-OP) 
and decouples the system into electromagnetically cou- 
pled OPs in the underdoped regime. 

4. Conclusions 

High quality single crystals of 

Bi2+a^Sr2-a:Ca2Cu30io+5 (Bi2223) with x = 0.2 were 
grown using a traveling solvent floating zone method 
in order to investigate dimensionality of superconduc- 
tivity in highly underdoped Bi2223 crystals. Obtained 
crystals were characterized by X-ray diffraction, DC 
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susceptibility and resistivity measurements, confirming 
Bi2223 to be the main phase. The highest quahty crystal 
was annealed under various oxygen partial pressures, 
adjusting its carrier density from slightly underdoped 
to highly underdoped. From the relationship between 
Tc and c-axis length, we could conclude that higher Bi 
substitution for Sr {x > 0.1) is successfully achieved. 
However, Bi2223 with x = 0.2 did not show a decrease 
in Tc,max7 compared with x = 0.1 Bi2223. This behavior 
confirms that of the inner layer limits the bulk in 
the optimally-doped region. The fluctuation diamagnetic 
component was extracted from the anisotropic normal 
state susceptibilities Xab{T) {H ± c) and Xc{T) {H \\ c). 
The temperature dependence of this component became 
strong with underdoping, suggesting a reduction of the 
superconducting dimensionality and/or an increase of 
vortex liquid domain. This behavior supports the view 
that inner Cu02 layer which is relatively underdoped 
compared to outer layers disrupts superconducting 
coherence among the three layers and change interlayer 
coupling situation more largely than that of Bi2212. 
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